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Nuncrical  Solution  of  Inverse  Problems 
Statement  of  Problems  Studied 

It  was  proposed  to  derive  efficient  3  dimensional  numerical  algorithms  for  in¬ 
version  of  model  equations:  in  geophysics,  for  underground  detection;  and  in  ultra 
sonic  tomography,  for  accurate  image  reconstruction  of  a  part  of  the  interior  of  a 
human  being. 

An  appropriate  model  for  both  of  the  above  problems  is  the  Helmholtz  equation 
model 


This  is  the  simplest  model  which  contains  all  of  the  difficulties  that  one  encounters, 
even  in  the  inversion  of  more  complicated  models,  such  as  the  Masowll  equations. 


The  problem  is  to  reconstruct  the  function  f  ,  based  on  measuring  u  on  the 
boundary  of  the  region. 

Stannary  of  Most  Inportant  Results 


they  are  also  optimally  efficient  based  on  mathematical  results  [1]  obtained 
uider  previous  ARO  sponsorship.  The  third  result  is  inportant  in  the  applications 
of  mathematics,  since  it  explains  how,  and  vhen  the  well-tacwn  rational  fmction 
algorithms  work.  _ . _ _ 


In  the  case  of  ultrasonic  tonography, 
u  denotes  the  sound  pressure,  c(r)  is 
the  "speed"  of  sotnd  in  the  body  B  ,  Cq 
is  the  speed  of  sotnd  in  the  meditxn 
(usually  a  liquid)  surrounding  B  ,  and 
at  denotes  the  frequency.  In  this  case 
u  -  u^  satisfies  the  relation 


(4)  *00  s  ^  log  u^(r)  |Id  -  |  /Tf 1  ds  +  0(k"°)  ,  k  -  - 

rs  P 

where  P  is  the  path  of  sotnd  fran  the  source  point  rg  to  the  detector  point  r^  . 

We  have  developed  an  effective  method  of  reconstructing  f  ,  based  on  any  model 
for  which  the  sotnd  pressure  u^  satisfies  (4) .  Our  method  predicts  *(»)  -  /p  /L+f  ds 
via  a  staple  rational  function  algorithm.  In  practice  u^CFg)  and  u^Cr^)  may 
be  measured  for  all  frequencies  (recall  k  -  u/cq  ,  u>  -  2ir  x  frequency)  in  (e.g.)  the 
range  1  megahertz  <  u/2ir  <  4  megahertz.  Our  experiments  have  shown  that  if  u^(rg) 
and  ^(Fj)  are  known  accurately,  then  we  can  accurately  reconstruct  f  via  well- 
known  x-ray  algorithms.  In  practice  one  has  noise,  and  there  is  error  in  the 
measurements.  We  have  thus  also  developed  a  simple  procedure  for  circumventing 
the  effects  of  these  errors.  We  axe  thus  able  to  distinguish  objects  having  a 
diameter  of  about  one  wavelength  and  separated  by  about  a  wavelength. 

A  publication  of  these  results  will  reach  your  office  before  Nov.  30,  1982. 
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52. 


The  author  has  also  developed  a  procedure  for  reconstructing  the  function  f(r) 
in  the  equation  (2)  in  the  region  {r  -  (x,y,z):  z  >  0}  ,  assuring  that  u^Cx.y.O) 
is  known  for  0  <  k  <_  «  .  This  latter  procedure  is  based  on  a  new  formula  for  f  , 
some  new  approximations,  and  sine  approximation  formulas  developed  by  the  author 
[1]  tzider  previous  ABO  support.  It  has  the  desirable  feature  that  it  requires  no 
matrix  inversion.  Rather,  the  3-dimensional  problem  "decouples"  and  requires  only 
the  numerical  approximation  of  one  and  two  dimensional  integrals.  A  paper  on  this 
procedure  will  reach  your  office  before  Dec.  20,  1982. 


53. 

Probably  the  author's  most  important  discovery  from  a  mathematical  standpoint 
is  in  the  area  of  rational  functions.  We  would  not  have  been  able  to  solve  the 
problem  SI  above  without  this  result.  We  first  describe  a  special  case  of  this 
result,  and  then  point  out  its  consequences. 

(a)  If  p  >  1  ,  q  «*p/(p-l)  ,  if  g  is  analytic  in  the  unit  disc  U  of 
the  cocplex  plane,  such  that 


(5) 


lhn  f  |G(rei6)  P  de) 

r+l~  J0 


where  G(z)  -  ,  then  by  taking 

1-z 

n  v2 

(6)  h  -  irOj)  ,  S(k,h).(x) 


sin[£  (x-kh)  ] 
£  (x-kh) 


(7) 


zk  -  (e^-D/Ce^+l) 


(8) 
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we  have 


(»)  IgW  -  ?  g(V  S(k,h).<K>Ol  <  CpN1/2  «p{-it(^)1/2> 


The  bound  on  die  right  hand  side  is  essentially  the  best  possible,  under  any  2N  +  1  - 
point  method  of  approximation. 

This  result  [1]  was  discovered  by  the  author  un^er 
previous  AR0  sponsorship.  While  this  result  was  reported  earlier,  we  mention  it 
here  again,  since  it  led  to,  and  is  related  to  the  results  (b),  (c)  vhich  follow. 

(b)  Under  the  same  conditions  in  (a)  above  one  sets 

2  N  z-z, 

(10)  p(z)  -  (l-z2)  n 

k--N  1  V 


then 

ai) 


N  p(x)  g(zj 

l8W  '  J-H 


|  <  r  N172  exp{-7r(^)1/2} 


That  is,  we  have  an  explicit  rational  approximation  of  g  for  which  the  error 
converges  at  the  same  optimal  rate  as  that  in  (9) .  Because  of  its  optimal  nature 
this  linear  rational  approodmation  procedure  is  very  nearly  as  good  as  the  best 
rational  approximation  of  g  ,  of  the  same  degree.  It  is  curious  that  the  (almost) 
optimal  interpolation  points  are  the  "sine"  points  (7) . 

(c)  Other  consequences  of  the  result  (b)  : 

(i)  It  is  easy  to  develop  a  Thiele  algorithm  which  yields  the 
rational  function  in  (11). 

(ii)  The  problem  of  poles  on  the  interval  of  approximation,  vhich  has 

plagued  users  of  rational  function  in  the  past,  has  been  eliminated, 
(ill)  We  have  described  a  space  of  functions  in  which  rational  approximation 
works.  We  can  thus  tell  a  priori,  when  e.g.  the  "Hxj-algorithm"  for 
prediction,  or  the  ODE  method  of  Bulirsh  end  Stoer  work. 
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(iv)  The  transformation 


is  a  conformal  map  of  the  unit  disc  onto  the  right  half  plane.  It 
also  transforms  the  rational  function  in  (11)  into  another  rational 
function  for  interpolation  over  (0,®)  ,  and  a  description  of  a 
space  of  functions  (e.g. ,  analyticity  and  boundedness  on  the  right - 
half  plane,  of  class  Lip^y^  on  (0,®))  for  which  this  works. 
The  corresponding  rational  function  approximation  is 


G(0 


N  i|/(OG<ek) 

J-n  te-y*‘«k> 


HO 


N 


where  ^  -  e1*  ,  and  the  error  is  OfN1^2  exp{-ir(^)  }]  ,  and  there 
we  have  also  assured  that  G(0)  *  G(®)  -  0  . 

(v)  Similary  £  ■  ew  yields  a  rational  function  of  ew  ;  it  describes  a 
space  of  function  for  thich  this  is  an  optimal  approximation,  and 


in  particular,  it  describes  a  space,  and  a  set  of  interpolation 


points  (w^  «  fch)  for  which  P.  Wynn's  "epsilon  algorithm"  -  which 
includes  the  Fade  method  -  always  works. 


We  hope  to  have  these  results  in  your  office  by  the  end  of  January,  1983. 

(vi)  Our  result  plays  the  same  role  for  rational  approximation  as  the 
Chebyshev  polynomials  play  for  approximation  by  polynomials.  If 
F  is  analytic  in  the  ellipse  E  with  foci  at  +  1  and  sun  of 
semi-axes  equal  to  R  (R  >  1)  then  Lagrange  interpolation  at  the 
nodes  of  the  Chebyshev  polynomial  T^  produces  very  nearly  as 
good  a  polynomial  of  degree  n  approximation  to  F  on  [-1,1] 
as  the  best  (uniform)  polynomial  approximation.  Since  the  best 


triform  polyncmlal  approximation  is  unduly  difficult  to  obtain,  one 
is  almost  always  satisfied  with  the  nearly  best  approximation  pro* 
vided  by  the  Chebyshev  polynomial. 
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